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Abstract  -  The  Moderate  Resolution  Imaging  Spectroradiometer 
(MODIS)  makes  observations  using  36  spectral  bands  with 
wavelengths  from  0.41  to  14.4pm  and  nadir  spatial  resolutions  of 
0.25km,  0.5km,  and  1km.  It  is  currently  operating  onboard  the 
NASA  Earth  Observing  System  (EOS)  Terra  and  Aqua  satellites, 
launched  in  December  1999  and  May  2002,  respectively.  The 
MODIS  Level  IB  (LIB)  program  converts  the  sensor’s  on-orbit 
responses  in  digital  numbers  to  radiometrically  calibrated  and 
geo-located  data  products  for  the  duration  of  each  mission.  Its 
primary  data  products  are  top  of  the  atmosphere  (TOA) 
reflectance  factors  for  the  sensor’s  reflective  solar  bands  (RSB) 
and  TOA  spectral  radiances  for  the  thermal  emissive  bands 
(TEB).  The  LIB  algorithms  perform  the  TEB  calibration  on  a 
scan-by-scan  basis  using  the  sensor’s  response  to  the  on-board 
blackbody  (BB)  and  other  parameters  which  are  stored  in  Look¬ 
up  Tables  (LUTs).  The  RSB  calibration  coefficients  are  processed 
offline  and  regularly  updated  through  LUTs.  In  this  paper  we 
provide  a  brief  description  of  the  MODIS  LIB  calibration 
algorithms  and  associated  LUTs  with  emphasis  on  their  recent 
improvements  and  updates  developed  for  the  MODIS  collection  5 
processing.  We  will  also  discuss  sensor  on-orbit  calibration  and 
performance  issues  that  are  critical  to  maintaining  LIB  data 
product  quality,  such  as  changes  in  the  sensor’s  response  versus 
scan-angle. 

I.  Introduction 

Two  Moderate  Resolution  Imaging  Spectroradiometer 
(MODIS)  instruments,  a  Proto-flight  Model  (PFM)  and  a  Flight 
Model  1  (FM1),  were  designed  and  built  by  Raytheon  /  Santa 
Barbara  Remote  Sensing  (SBRS),  Goleta,  California,  for  the 
NASA’s  Earth  Observing  System  (EOS)  missions  that  promote 
the  continuity  of  essential  science  and  environmental  data 
records  [1-5].  The  MODIS  PFM  onboard  the  Terra  spacecraft 
has  been  in  operation  for  more  than  5  years  since  its  launch  in 
December  1999  and  the  FM1  for  nearly  3  years  onboard  the 
Aqua  spacecraft  launched  in  May  2002.  Each  instrument  makes 
continuous  observations  over  a  wide  field-of-view  (FOV)  in  36 


spectral  bands  with  wavelengths  from  0.41  to  14.4pm  and  at 
nadir  spatial  resolutions  of  0.25km,  0.5km,  and  1km.  Table  1 
lists  key  sensor  specified  parameters  and  characteristics. 


Table  1  MODIS  Specified  Parameters  and  Characteristics 


Orbit 

705km  (altitude) 

Sun-synchronous,  near-polar,  circular 

10:30AM  descending  node  (Terra) 

1 :30PM  ascending  node  (Aqua) 

Spectral  Range 

20  reflective  solar  bands  (1-19  and  26)  from 

0.41  to  2.2p, 

16  thermal  emissive  bands  (20-25  and  27-36) 
from  3.7  to  14.4p 

Spatial  Resolution 
(nadir) 

0.25km  (bands  1-2,  40  detectors  along-track  and 

4  sub-samples  each  cross-track) 

0.5km  (bands  3-7,  20  detectors  along-track  and 

2  sub-samples  each  cross-track) 

1km  (bands  8-36,  10  detectors  along-track) 

Swath 

10  (1km  resolution)  pixels  along-track  by  1354 
(1km  resolution)  pixels  cross-track 

Scan  Rate 

20.3  rpm  ( 1 .477s  /  scan) 

Scan  Range 

±55°  from  nadir  or  angle  of  incident  (AO I)  to 
the  scam  mirror  from  10.5° to  65.5° 

Digitization 

12  bits  ADC 

Data  Rate 

1 1Mbps  daytime 

3.3Mbps  nighttime 

On-board 
Calibrators  (OBCs) 

Solar  diffuser  (SD) 

Solar  diffuser  stability  monitor  (SDSM) 
Blackbody  (BB ) 

Spectro-radiometric  calibration  assembly 
(SRCA) 

Focal  Plane 
Assemblies  (FPAs) 

Visible  (VIS) 

Near  infrared  (NIR) 

Short-  and  mid-wave  infrared  (SMIR) 

Long- wave  infrared  (LWIR) 

Spectral  Bands  and 
Detectors 

36  bands 

490  detectors  (including  time-delay  and 
integration  detectors  used  for  the  high/low  gain 
bands  13  and  14) 

The  MODIS  Level  IB  (LIB)  algorithms  convert  the 
sensor’s  responses  in  digital  numbers  (DN)  to  top  of  the 
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atmosphere  (TOA)  reflectance  factors  for  the  reflective  solar 
bands  (RSB)  and  TOA  spectral  radiances  for  the  thermal 
emissive  bands  (TEB).  They  are  designed  based  on  instrument 
performance  determined  from  comprehensive  pre-launch  and 
on-orbit  calibration  and  characterization  [6-9].  Two  separate 
versions  of  the  LIB  software  and  corresponding  Look-up 
Tables  (LUTs),  one  for  Terra  MOD1S  and  the  other  for  Aqua 
MOD1S,  are  maintained  to  accommodate  differences  between 
the  two  sensors  and  to  track  their  on-orbit  changes 
independently.  The  NASA  MODIS  Characterization  Support 
Team  (MCST)  is  responsible  for  the  tasks  of  performing  sensor 
calibration  and  characterization,  maintaining  and  improving  the 
LIB  algorithms,  and  deriving  and  updating  the  calibration 
parameters  in  the  LIB  LUTs.  In  this  paper  we  provide  a  short 
description  of  the  MODIS  LIB  calibration  algorithms  and 
associated  LUTs  and  highlight  the  latest  improvements  and 
updates  implemented  for  the  MODIS  collection  5  processing. 

II.  Overview  of  MODIS  LIB  Calibration  Algorithms 

The  MODIS  LIB  calibration  is  performed  for  each  band, 
detector,  sub-sample,  and  mirror  side  (BDSM).  The  algorithm 
consists  of  two  major  modules:  one  for  the  reflective  solar 
bands  (RSB)  and  the  other  for  the  thermal  emissive  bands 
(TEB).  Details  of  the  MODIS  TEB  and  RSB  calibration 
algorithms  are  provided  in  papers  given  by  Xiong,  et  al  [10- 
11].  The  LIB  input  includes  the  L1A  data  (raw  counts),  geo¬ 
location  data,  and  LIB  LUT  files  for  the  TEB,  RSB  and  quality 
assurance  (QA).  A  simplified  LIB  data  processing  flow 
diagram  is  shown  in  Figure  1. 
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Figure  1  MODIS  Level  IB  (LIB)  Flow  Diagram 


The  thermal  emissive  bands  (20-25  and  27-36)  are 
calibrated  by  an  on-board  blackbody  on  a  scan-by-scan  basis 
using  a  quadratic  algorithm  that  corrects  for  instrument 
background  (from  the  sensor’s  space  view),  thermal  emissions 
from  the  scan  mirror  and  scan  cavity,  and  response  versus  scan 
angle  (RVS).  An  optical  leak  correction  is  applied  to  bands  32- 
36  (Terra  only).  The  retrieved  Earth  view  (EV)  radiance  is 
converted  to  a  scaled  integer  (SI). 


The  reflective  solar  bands  (1-19  and  26)  are  calibrated 
using  a  solar  diffuser  (SD)  and  solar  diffuser  stability  monitor 
(SDSM)  system.  The  SD/SDSM  calibration  is  typically 
performed  bi-weekly  (weekly  for  the  first  year)  and  analyzed 
off  line  by  the  MCST.  The  RSB  calibration  uses  a  linear 
algorithm  with  calibration  coefficients  updated  through  the 
LUTs.  The  algorithm  corrects  for  instrument  background, 
temperature,  and  response  versus  scan  angle  (RVS).  A  thermal 
leak  correction  is  applied  to  the  SWIR  bands  (5-7  and  26).  The 
retrieved  EV  reflectance  factor  is  converted  to  a  scaled  integer 
(SI).  The  RSB  calibration  is  executed  after  the  TEB  calibration. 

The  LIB  output  consists  of  calibrated  radiances  and 
reflectance  factors  and  is  organized  in  three  Hierarchical  Data 
Format  (HDF)  files  corresponding  to  MODIS’  three  spatial 
resolutions  and  their  associated  metadata  files.  The  following 
expressions  are  used  to  convert  the  SI  to  the  scene  TOA 
radiances  or  reflectance  factors. 

For  TEB  radiances: 

Radiance  =  radiance  _scales  *  (SI  —  radiance  _offsets) 

For  RSB  reflectance  factors: 

Reflectance  =  reflectance  _scales*  (SI  -reflectance  ^offsets) 

For  RSB  radiances: 

_  ,.  _  Solar  Irradiance  over  pi 

Radiance  =  Reflectance - = - = - =f— 

( Earth  -  Sun  _  Distance ) 

The  band  dependent  coefficients  of  radiance_scales, 
radiance _offsets,  reflectance _scales,  and  reflectance  ojfsets 
are  stored  as  scientific  data  set  (SDS)  attributes  within  the 
Level  IB  output  files.  The  “Earth-Sun  Distance”  and 

“ Solar _lrradiance_over_pr  are  provided  as  global  attributes. 
The  LIB  also  produces  a  separate  file  containing  on-board 
calibration  (OBC)  data  and  key  instrument  telemetry.  The 
OBC  file  is  primarily  used  for  off-line  calibration  and 
characterization  purposes. 

III.  MODIS  LIB  Look-Up  Tables 

The  MODIS  Level  IB  (LIB)  code  Look-up  Tables  (LUTs) 
are  organized  into  three  groups:  one  for  the  TEB  calibration, 
another  for  the  RSB  calibration,  and  a  third  one  for  quality 
assurance  (QA).  Each  individual  LUT  is  implemented  in  the 
appropriate  HDF  file  as  either  a  global  attribute  or  a  Scientific 
Data  Set  (SDS)  [12],  Currently  there  are  98  LUTs  (40  TEB,  34 
RSB,  and  24  QA)  for  the  Terra  MODIS  LIB  code  and  101 
LUTs  (43  TEB,  34  RSB,  and  24  QA)  for  the  Aqua  MODIS 
LIB  code.  Aqua  MODIS  has  3  extra  TEB  LUTs  compared 
with  Terra  MODIS.  These  extra  LUTs  are  used  to  provide  pre¬ 
determined  linear  gain  coefficients  for  bands  33,  35,  and  36  at 
high  BB  temperatures  during  BB  warn-up  process  when  their 
responses  saturate  (a  feature  that  only  appears  in  the  Aqua 
MODIS). 


Some  of  the  LUTs  included  in  the  HDF  files  are  not  used 
for  calibration,  such  as  those  used  for  serial  numbers,  metadata 
items,  and  code  version  numbers.  Some  of  the  LUTs  are 
constant  and  some  may  depend  on  the  data  collection  time  in 
either  a  “step  function’’  or  “piecewise  linear”  style.  Among 
calibration  LUTs,  some  were  derived  from  pre-launch 
calibration  and  characterization,  such  as  the  relative  spectral 
response  (RSR),  the  temperature  correction  coefficients,  and 
corresponding  reference  temperatures.  Some  calibration  LUTs 
are  derived  from  on-orbit  measurements  and  updated  as 
needed.  Examples  include  the  RSB  response  versus  scan  angle 
(RVS)  and  SD  calibration  coefficients  for  the  reflectance  factor 
retrieval.  The  MCST  at  NASA/GSFC  is  responsible  for 
determining  the  contents  of  all  LIB  LUTs  and  for  maintaining 
and  updating  (adding,  deleting,  or  changing)  them. 

IV.  Status  OF  MODIS  LIB  and  LUTs 

Terra  MODIS  has  been  in  operation  for  more  than  five 
years  since  its  launch  in  December  1999.  Thus  far  there  have 
been  four  major  configurations  used  for  the  instrument 
operation.  As  shown  in  Table  2  that  there  are  12  versions  of  the 
LIB  code  used  in  data  production.  Each  change  was  made  to 
improve  the  LIB  algorithms  or  to  mitigate  code  deficiencies. 
Many  versions  of  the  LUTs  have  been  updated  either  to  the 
track  sensor’s  response  changes  over  time  or  to  accommodate 
changes  in  operational  configuration. 


Table  2  Terra  MODIS  LIB  code  and  LUT  versions 


Code  Version 

Production  Time 

LUTs  Version 

V2.3.2 

03/17/2000  (077  2000)  00:00 

3 

V2.4.2 

06/19/2000(171  2000)  00:00 

0 

V2.4.3 

08/18/2000(231  2000)  14:00 

1 

V2.4.4 

10/13/2000  (287  2000)  19:55 

0-2 

V2.5.4 

11/23/2000  (328  2000)  15:55 

0 

V2.5.5 

02/13/2001  (044  2001)  13:55 

1,2 

V3.0.0 

05/24/2001  (144  2001)  00:00 

2-7 

V3.0.1 

02/25/2002  (056  2002)  00:00 

0,  1 

V4.1.2 

01/30/2003  (030  2003)  01:55 

1-11 

V4.2.0 

08/22/2003  (234  2003)  02:00 

3-9 

V4.3.0 

12/22/2003  (356  2003)  22:35 

1-9,  11-25 

V5.0.6 

03/07/2005  (066  2005)  23:55 

1,2 

Since  its  launch  in  May  2002,  Aqua  MODIS  has  been  in 
operation  for  nearly  three  years.  Unlike  its  Terra  predecessor, 
the  Aqua  MODIS  has  been  operated  using  the  same  B-side 
electronic  configuration  for  the  entire  mission.  Table  3  is  a  list 
of  its  LIB  code  and  LUT  versions. 


Table  3  Aqua  MODIS  LIB  code  and  LUT  versions 


Code  Version 

Production  Time 

LUTs  Version 

V3.1.0 

06/07/2002  (158  2002)  18:10 

0-3 

V4.1.1 

10/31/2002  (304  2002)  00:15 

0,1 

V4.1.3 

01/22/2003  (022  2003)  09:55 

0-11 

V4.2.1 

08/21/2003  (233  2003)  12:00 

4-8 

V4.3.1 

06/02/2002  (183  2002)  15:25 

1-16 

V5.0.7 

Scheduled  for  May  2005 

1 

The  latest  MODIS  data  collection  will  be  produced  using 
version  5  of  the  LIB  code.  The  Terra  collection  5  (V5.0.6)  was 
started  in  March  2005  and  the  Aqua  collection  5  will  start  in 
May  2005.  In  addition  to  the  NASA  Goddard  Distributed 
Active  Archive  Center  (GDAAC),  selected  LIB  products  have 
also  been  produced  at  the  MODIS  Adaptive  Processing  System 
(MODAPS)  and  at  the  MODIS  Ocean  Color  Data  Processing 
System. 

Only  a  few  minor  code  changes  were  made  in  LIB  version 
5,  including  improving  code  portability,  enabling  B21 
calibration  with  mirror  side  dependency  (for  reducing  image 
striping),  and  adding  new  LUTs  for  SWIR  correction  with 
flexible  detector  dependency  in  order  to  avoid  using  noisy 
detectors  in  the  correction  algorithm.  Overall  improvements  are 
primarily  coming  from  better  calibration  and  LUTs  that  are 
developed  using  long-term  consistent  calibration  data. 

In  version  5  the  Terra  TEB  will  use  the  response  versus 
scan-angle  (RVS)  derived  from  on-orbit  spacecraft  pitch 
maneuvers  made  on  March  26  and  April  14,  2003  for  the  entire 
mission,  replacing  the  relative  RVS  determined  from  scanning 
the  closed  nadir  aperture  door.  For  the  reflective  solar  bands  of 
both  Terra  and  Aqua  MODIS,  on-orbit  RVS  changes  of  mirror 
side  one  are  updated  using  each  sensor’s  response  trending 
from  SD  calibrations  and  lunar  observations.  The  RSB  RVS 
changes  of  mirror  side  two  are  derived  relative  to  mirror  side 
one  using  the  trending  of  mirror  side  response  ratios  at  various 
angles  of  incidence  (AOI).  Several  QA  LUTs  have  been 
updated  based  on  the  detectors’  on-orbit  performance. 
Complete  histories  of  LIB  code  versions  and  LUT  updates  for 
both  Terra  and  Aqua  MODIS  are  posted  on  the  MCST  web 
page  (http://www.mcst.ssai.biz/mcstweb/LlB/  product. html). 

Other  documents  that  describe  the  algorithms,  function  and 
applications  of  the  LIB  code  and  LUTs  include  “MODIS  Level 
IB  Algorithm  Theoretical  Basis  Document”,  “MODIS  Level 
IB  Product  User’s  Guide”,  “MODIS  LUT  Information  Guide”, 
“MODIS  Level  IB  In-Granule  Calibration  Code  High  Level 
Design”,  and  “MODIS  Level  IB  Products  Data  Dictionary”. 
These  are  also  provided  at  the  MCST  web  page  (home  page 
http://www.mcst.ssai.biz/mcstweb/). 

V.  Summary 

Since  their  launches  in  December  1999  and  May  2002, 
Terra  and  Aqua  MODIS  together  have  produced  nearly  eight 
years  of  data  products  for  the  science  community’s  studies  of 
the  Earth  system  and  its  land,  oceans,  and  atmosphere. 
Working  together  with  the  MODIS  Science  Team  members 
and  sensor  vendor  representatives,  the  NASA  MODIS 
Characterization  Support  Team  (MCST)  is  responsible  for 
providing  instrument  on-orbit  calibration  and  characterization, 
maintaining  and  improving  the  LIB  algorithms,  and  updating 
the  calibration  parameters  in  the  LIB  LUTs.  Currently  two 
separate  versions  of  the  LIB  software  and  LUTs  have  been 
maintained  and  updated  for  the  data  processing,  one  for  the 
Terra  and  the  other  for  the  Aqua  MODIS.  This  paper  has 
presented  the  status  of  the  MODIS  LIB  algorithms  and 
associated  LUTs,  including  the  latest  improvements  and 
updates  made  for  the  MODIS  collection  5  processing. 
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